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Rolling Balls and Octonions 



A. A. Agrachev* 



Abstract 

In this semi-expository paper we disclose hidden symmetries of a 
classical nonholonomic kinematic model and try to explain geometric 
meaning of basic invariants of vector distributions. 

1 Introduction 

The paper is written to the 70th anniversary of Vladimir Igorevich Arnold. 
This is just a small mathematical souvenir but I hope that Vladimir Igorevich 
will get some pleasure looking over it. The content of the paper is well 
described by the cryptogram below. The figure represents the roots system 
of the exceptional group G2 (the automorphisms group of octonions) and two 
touching each other circles whose diameters are in ratio 3:1. 




Our starting point is a classical nonholonomic kinematic system that is 
rather important in robotics: a rigid body rolling over a surface without 
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slipping or twisting. The surface is supposed to be the surface of another 
rigid body so that the situation is, in fact, symmetric: one body is rolling 
over another one. We also assume that the surfaces of the bodies are smooth 
and cannot touch each other in more than one point. 

This system has a 5-dimensional configuration space: coordinates of the 
points on each surface where the surfaces touch each other give 4 parameters; 
fifth parameter measures mutual orientation of the bodies at the touching 
point. 

Now let us imagine one of the bodies immovable and another rolling. 
Given an initial configuration, it is rather clear that one can roll the body 
in a unique way along any curve on the immovable surface starting from 
the initial touching point. In other words, given an initial configuration, 
admissible motions are parameterized by the curves on the two-dimensional 
surface with a fixed initial point. On the other hand, it is not hard to prove 
(see [2j Ch. 24]) that admissible motions allow to reach any configuration from 
any other as soon as two bodies are not congruent. We are thus in a typical 
nonholonomic situation since the configuration space is 5-dimensional. 

Now turn to Mathematics. Admissible velocities (i.e. velocities of admis- 
sible motions) form a rank 2 vector distribution on the configuration space. 
This distribution is involutive if and only if our bodies are the balls of equal 
radii. More precisely, let vector fields /, g form a local basis of our distribu- 
tion, [/, g] be the commutator (Lie bracket) of the fields, and q a point of 
the configuration space. It happens that vectors f(q),g(q), [f, g](q) are lin- 
early independent if and only if the curvatures of two surfaces are not equal 
at their touching points corresponding to the configuration q. Moreover, if 
these curvatures are not equal, then: 

/(<?) A g(g) A [f,g](q) A [f, [f,g]](q) A [g, [g,f]](q) ? 0. (1) 

In other words, the basic fields, their first and second order Lie brackets are 
all linearly independent. 

The germs of rank 2 distributions in M 5 with property (1) were first 
studied by Elie Cartan in his famous paper [6]. A brief count of parameters 
demonstrates that 5 is the lowest dimension where the classification of generic 
germs of distributions must have functional invariants. This happens for rank 
2 and rank 3 distributions in M 5 ; moreover, the classifications for both values 
of the rank are essentially equivalent and it is sufficient to study the case of 
rank 2. 
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Cartan found a fundamental tensor invariant that is a degree 4 symmetric 
form on the distribution: the desired functional invariant is just the cross- 
ratio of the roots of this form. There is exactly one equivalence class for which 
the Cartan's form is identical zero. We call the germs from this class flat; a 
germ of the distribution is flat if and only if it admits a basis generating the 
5-dimensional nilpotent Lie algebra. Cartan showed that symmetry group of 
a flat distribution is the split form of the 14-dimensional exceptional simple 
Lie group Gi- In all other cases, dimension of the symmetry group does not 
exceed 7. 

How all that is related to the rolling bodies? It is rather obvious that the 
group of symmetries of the "no slipping no twisting distribution" is trivial 
for generic bodies. This group acts transitively on the configuration space 
if and only if the surfaces of both bodies have constant curvatures. Let us 
stay with the case of constant nonnegative curvatures so that the surfaces are 
spheres (one of them may be a plane, i.e. the sphere of the infinite radius). 

Natural symmetries are isometries of the spheres; they form a 6-dimensi- 
onal Lie group. It was Robert Bryant who claimed first that no slipping no 
twisting distribution is flat in the case of spheres whose radii are in ratio 3:1. 
He insisted that just followed the Cartan's method and never published this 
fact in a paper or a preprint. It also remained absolutely unclear what are 
hidden additional symmetries (dimension of the symmetry group jumps from 
6 to 14 when ratio of the radii gets 3:1!). Unfortunately, Cartan's method 
does not give much in this regard: the construction of the fundamental in- 
variant is based on the involved reduction-prolongation procedures in the jet 
spaces, and any connection with the original geometric problem is lost quite 
far from the end of the way. 

The goal of this note is to finally untwine this puzzle. In Section 2 we 
give a simple twistor model for the configuration space of the rolling balls 
problem and for the no slipping no twisting distribution; the role of the group 
G2 is not yet clear from this model. 

Split form of G2 is the group of automorphisms of split- octonions, a hy- 
perbolic version of usual octonions, where positive definite quadratic form 
(square of the norm of octonions) is substituted by a non-degenerate sign- 
indefinite quadratic form. Split-octonions have nontrivial divisors of zero 
(zero locus of the sign-indefinite form). A simple quadratic transformation 
demonstrates that no slipping no twisting distribution is equivalent to the 
"divisors of zero distribution" of split-octonions if the radii of the balls are 
in ratio 3:1. This is the subject of Section 3. 
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In Section 4 we outline a recently developed variational way to construct 
differential invariants of vector distributions in order to put the rolling bodies 
model in a broader framework and to explain the meaning of basic invariants. 

Acknowledgements. I am grateful to Robert Bryant, Aroldo Kaplan and Igor 
Zelenko for very interesting stimulating discussions. 

2 Twistor model 

We study admissible motions of two balls B r and B R rolling one over another 
without slipping or twisting. Here r and R are the radii of the balls. Instant 
configuration of the system of two balls is determined by an orientation 
preserving isometry of the tangent planes to the spheres S r = dB r and Sr = 
8Br at the points where the balls touch each other. In other words, the state 
space of our kinematic system is: 

w / rp q , qi e S r , q 2 e S R , u. is an \ 

MR, r = < H . l qi b r — > l q2 b R . . >. 

[ isometry of oriented planes J 

It is easy to see that Mr^ is a smooth 5-dimensional manifold. Motions 
of the system are families of isometries /i(t) : T qi ^S r — > T Q2 ^Sr, tel. The 
no slipping condition reads: 

fx(t)( qi (t))=q 2 (t). 

The no twisting condition requires that /i(t) transforms parallel along q± (t) 
vector fields in the parallel along (^(O vector fields. 

These two conditions define a rank 2 vector distribution D R,r on Mr^. 
We have: D R,r = |J D R ' r , where D R,r is a two-dimensional subspace of 

T^Mr^; admissible motions of the two balls system are exactly integral curves 
of the distribution D R,r . Given an initial configuration, the ball B r can be 
rolled in a unique way along any smooth curve on Sr and the same is true if 
we transpose r and R. In the formal geometric language this observation just 
means that the subspace D R,r C T^Mr^, where u, : T qi — > T q2 , is projected 
one-to-one onto T qi SR and T q2 S r , V/x G Mr^. 

In what follows, we treat S r and Sr as unit spheres in R 3 with re-scaled 
metrics: 

M R , r = {/i : qi -> q^ \ % G M 3 , |%| = 1, i = 1,2, ^2|a*(«)| = ?>|, G 
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Let p — j, the homothety i p : p i— > p/x, p G Mr )7 ., transforms M^ r in M^i. 
We set D p = L p *D R,r . The distribution D p on M^i is determined by the 
"re-scaled no slipping condition" 

= pg 2 (i) 

and the no twisting condition; the last one remains unchanged. 

From now on we will deal with fixed space M 1;1 endowed with the family 
of distributions D p instead of the family of pairs (M RjT , D R ' r ) . In order 
to explicitly describe the distributions D p we use a classical quaternion's 
parameterization of spherical bundle p : S — > S 2 , to the unit sphere S 2 Ci 3 , 
where 

S = {(q, v) G R 3 x R 3 I \q\ = \v\ = 1, (g, u) = 0}, p(g, v) = q. 

Let us recall this parameterization. We identify R 3 with the space of imagi- 
nary quaternions: 

M 3 = {ai + f3j + 7/c | 7 G 1R} C H. 

Let S 3 = {to G H |w| = 1} be the group of unitary quaternions; then 
f) : S 3 — > S* 2 , h(w) = tuuu is the classical Hopf bundle, while the mapping 

\& : S* 3 — > <S, ^(iu) = (wiw,wjw) 

is a double covering. Moreover, the diagram 

S 3 «S 
fj p 
S 2 S 2 

is commutative, hence \1/ is a fiber-wise mapping of the bundle t) : S 3 S 2 

onto the bundle p : 5 — ► 5 which induces a double covering of the fibers. 
The fibers of the bundle f) : S 3 — > S* 2 are the residue classes {e* e u> | G 
]Rmod27r} of the one-parametric subgroup generated by i. 
The normal to the fibers distribution 

span {jw, kw} C T W S 3 , w G S 3 
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is a connection on the principal bundle f) : S 3 — > S* 2 . It is easy to see that 
\1/ transforms this connection in the Levi Civita connection on the bundle 
p : S — > S* 2 , which defines standard parallel transport on S* 2 . 

We are now ready to give a quaternion's model of the rolling balls config- 
uration space Mi ; i (more precisely, of a double covering to M^i) and of the 
no slipping no twisting distributions -D p . For any Wi,w 2 G S* 3 there exists a 
unique orientation preserving isometry of the fiber 

= {e if V | G Mmod27r} 

on the fiber t)~ 1 {t){w 2 )) that sends iUi to w 2 . This isometry sends e* 9 u>i to 

e w 2 . Moreover, pairs (iui, ■u; 2 ) and (u^, u> 2 ) define one and the same isometry 

if and only if w[ = e %e w±, w 2 = e ld w\ (with one and the same 9'). Hence 

the coset space of S 3 x S* 3 by the action (wi,w 2 ) 1— > (e l6 wi, e t9 w 2 ) of the 

one-parametric group {e l6 \ 9 G Mmod27r} is a double covering of M 1:L . We 

use symbol M for this coset space and tt : S 3 x S* 3 — > M for the canonical 

projection. By D p = |J D£ we denote the rank 2 distribution on M that is 

xeD 

the pullback of D p by the double covering M — > M^i. Then 

D £( TO1 ,™ 2 ) = span {(jw u pjw 2 ), (kw u pkw 2 )} Vwi, w 2 G S 3 . 

Let us now treat the quaternionic space H 2 = {(iu^u^) | G H} as C 4 , 
where u>i = z\ + z 2 j, w 2 = z 3 + z 4 j, zjGC, I — 1, . . . , 4. We see that M is 
nothing else but a complex projective conic, 

M = { Zl :z 2 :z 3 :z 4 \ \ Zl \ 2 + |^ 2 | 2 = \z 3 \ 2 + | ^ 4 1 2 > C CP 3 . 

This conic is often called "the space of isotropic twistors" . Moreover, 

K(w uW2 ) =n*Cj(w 1 , P w 2 ). 

3 Split-Octonions 

We also treat H 2 as the algebra O 
where 

(a + £b)(c + £d) = 



6 



= {wi + £w 2 \ Wi EM} of split-octonions, 
= (ac + db) +£(ad + cb). (2) 



Let x = w 1 + £w 2 , x = W\ — £w 2 and Q(x) = xx = \wi\ 2 — |w 2 | 2 - Then 
Q(xy) = Q(x)Q(y) and x~ x = as soon as Q(x) ^ 0; the cone Q _1 (0) 
consists of the divisors of zero. 

The automorphisms group of the algebra O is the split-form of the excep- 
tional Lie group G 2 . These automorphisms preserve quadratic form Q and 
hence its polarization Q(x,y) = \{Q{x + y) — Q(x — y)). In particular, these 
automorphisms preserve the subspace R 7 = {x G O Q(l,x) = 0} and the 
conic 

K = {x e R 7 1 Q(x) = 0} = {x e 6 | xx = 0}. 

Moreover, the automorphisms group of O acts transitively on the "spheriza- 
tion" K = {R + x \ xeK\0} = S 2 xS 3 ofthe cone K. 

To any x E K \ there associated a 3-dimensional subspace of divisors 
of 

A x = {y E R 7 | xy = 0}; 

the "spherization" turns A x into a 2-dimensional subspace A x C T X K, where 
x = M. + x. Obviously, the automorphisms group of O preserves vector distri- 
bution A = {A x } xeK - 

Proposition 1 The mapping $ : (wi + £w 2 ) i— > [w^ l iwi +£(wi 1 w 2 )) induces 
the diffeomorphism of M onto K. Moreover, the differential of this diff'eo- 
morphism transforms the "no slipping, no twisting" distribution D 3 in the 
"divisors of zero" distribution A. 

Proof. Let $ : M — > K be the mapping induced by $. We give an explicit 
formula for take v 1 G S 2 , v-y = f)(wi), then 

+ iv 2 ) = 7r(wi + ^(wiu 2 )), Vw 2 G 5 3 . 

Now we have to prove that <&(x) (D x &y) = for any x = ui\ + £w 2 , y = 
zjw\ + 3£(zjw 2 ) such that |iui| = |u> 2 |, z G C. It is sufficient to make 
calculation in the case |u>i| = \w 2 \ = 1. We have: 

D Wl+iW2 ® {zjwi + £{zjw 2 )) = 2w 1 zkw 1 + 2£(w 1 zjw 2 ). 

The desired result now follows from the multiplication rule (2). 

Let us give an explicit parameterization of the distribution A on K. First 
we parameterize K itself: 

K = {«! + £v 2 | vi G M 3 , v 2 G M, |ui| = |u 2 | = !}■ 
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Then A Vl+ £ V2 = {viu+£(uv 2 ) \ u, {viu) G M 3 }; this is a simple corollary of the 
multiplication rule (2). Let v G M 3 , \v\ = 1; the mapping w i— > w + iwi; maps 
H onto the subspace {u \ u, (vu) G M 3 }. Now we substitute u by w + v\WV\ 
in the above description of A Vl+ £ V2 and obtain the final parameterization: 

A Vl+£v2 = {[vi,w] + £((w + viwvi)v 2 ) | w G H}. 

4 Jacobi curves 

In this section we briefly describe the variational approach to differential 
invariants of vector distributions (see [H [3l [H]) in order to put the rolling 
balls model in a wider perspective. This approach is based on the contempo- 
rary optimal control techniques and suggests an alternative to the classical 
equivalence method (see [51 El EI]) of Elie Cartan. 

Rank k vector distribution A on the n-dimensional smooth manifold M 
is just a smooth vector subbundle of the tangent bundle TM: 

A=(jA 9 , A q C T q M, dimA g = A;. 

Distributions A and A' are called locally equivalent at q G M if there exists 
a neighborhood O qo C M of q and a diffeomorphism $ : O qo O qo such 
that $,A 9 = A' Hq) , Mq G O qo . 

A local basis of A is a fc-tuple of smooth vector fields fi, ■ ■ ■ , fk G Vec M 
such that 

A q = spanifciq), . . . , f k (q)}, q e O qo . 
Given a local basis, one may compute the flag of the distribution: 

A l q = span {(ad4 • • • ad/-J io )(g) | < j < 1} , 1 = 1,2,..., 

where ad/ g d = [f,g] is the Lie bracket. 

It is easy to see that the subspaces A' do not depend on the local basis. 
We set A' = [J A^, a growing sequence of subsets in TM. This sequence 

stabilizes as soon as A' +1 = A'. The distribution is involutive if and only if 
A 2 = A, it is completely nonholonomic (our subject) if A' = TM for suffi- 
ciently big I. Generic distributions are, of course, completely nonholonomic. 
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Integral curves of the distribution are often called horizontal paths. It is 
convenient to consider all paths of class H 1 , not only smooth ones. We thus 
have a Hilbert manifold fi A of horizontal paths: 

= {7 e H\[0, 1]; M) I 7(i) G A 7(t) ,for a.e. t G [0, 1]}. 

Now consider boundary mappings 

d t :tt A ^ M x M 

defined by the formula d t : (7) = (7(0), j(t)). It is easy to show that d t are 
smooth mappings. 

Critical points of the mapping d\ are called singular curves of A. Any 
singular curve is automatically a critical point of d t , Vi G [0,1]. Moreover, 
any singular curve possesses a singular extremal, i. e. an if 1 -curve A : [0, 1] — > 
T*M in the cotangent bundle to M such that 

A(t)GT 7 * w M\{0}, (A(t),-A(0)) J D 7 ^ = 0, Vie [0,1]. 

We set: 

A^ = {v G r 9 *M I ^,A 9 ) = 0,1/^0}, A ± = |J 

ijGM 

Obviously, A -1 is a smooth (n + /c)-dimensional submanifold of T*M (the 
annihilator of A). 

Let a be the canonical symplectic structure on T*M. Pontryagin Maxi- 
mum Principle implies that a curve A in T*M is a singular extremal if and 
only if it is a characteristics of the form <r| A ±; in other words, 

X(t) G A -1 , A(*) G ker(cr| A± ) , < t < 1. 

All singular extremals are contained in the characteristic variety 

C A = {z G A 1 - I kera 2 | A± ^ 0} . 

We have: Ca = A 2± if = 2; C A = A- 1 if k is odd; typically, C A is a codim 
1 submanifold of A 1 - if is even. 

A complete description of singular extremals is a hard task; to simplify 
the job we focus only on the regular part of the characteristic variety. We 
set: 

C A = {z£ Ca I dim kera z \ A± < 2, dim kero^x n T z Ca = 1} • 
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If k = 2, then C% = A 2± \ A 3± . 

C° A is a smooth submanifold of A- 1 ; it is foliated by singular extremals 
and by the fibers T*M n C° A . 

The motion along singular extremals defines a local flow on C\\ typically 
this flow is not fiber-wise, i. e. it does not transform fibers into fibers (see 
the figure above). 

Given z G C\ and an appropriate small neighborhood C z of z in we 
consider the canonical projection 

F : C z — ► C°/{singular extremals foliation} 

of C z on the space of contained in C z singular extremals. 

Assume that A is a singular extremal through z and it is associated to 
a singular curve 7, i.e. A(0) = z, X(t) G T*^M. Consider a family of 
subspaces 

J° x (t)=T x F(T; {t) MnC° z ) 

of the space 

T A C°/{singular extremals foliation} T Z C°JT Z X. 

Then t 1— > J°(£) is a smooth curve in the correspondent Grassmann mani- 
fold. Geometry of the curves t 1— > reflects the dynamics of the fibers 
along singular extremals and contains the fundamental information about 
distribution A. 

In what follows we assume that k = 2, n > 5 and A^, A 3 have maximal 
possible dimensions, i.e. dimA^ = 3, dim A 3 = 5. 
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i. First we consider the case n = 5 that is the case studied by Cartan (see 
Introduction). Let z G T*M and vr : T Z (T*M) -> T g M be the differential at 
z of the projection T*M -> M; then vr(J^(t)) C z x C T 9 M. Moreover, T g 7 C 
7r(J°(t)) and 1 1 — >■ tt (J\(t)) is a curve in the projective plane P (z ± /T q j). 

Proposition 2 (see [4]) A rank 2 distribution A on the 5- dimensional man- 
ifold is flat if and only if the curve 7r(J°(-)) is a quadric, for any singular 
extremal A. . 

In general, 7r ( is not a quadric. Let K z (q) C z -1 be the best approx- 
imating quadric for the curve 7r(J°(-)) near the point corresponding to the 
zero value of parameter t (the osculating quadric of classical projective geom- 
etry); then K z (q) is zero locus of a signature (2, 1) quadratic form on z /T q "f. 
We can, of course, treat K z (q) as zero locus of a degenerate quadratic form 
on z L . Finally, JC(q) = \J K z (q) is zero locus of a (3, 2) quadratic form on 

T 9 M (see 0] for details). 

The family of quadratic cones JC(q), q G M is an intrinsically "raised" 
from A conformal structure on M and A g C K(q). This conformal structure 
was first found by Nurowski [TT] who used the Cartan's equivalence method. 

Remark. Nurowski's conformal structure has a particularly simple descrip- 
tion for the divisors of zero distribution A of Sec. 3. Namely, /C(x) = 
n T x K , x G K, in this case. 

ii. From now on n is any integer greater or equal to 5. Let z G C^, A be the 
singular extremal through z and 7 the correspondent singular curve. We set: 

J x (t) = D\F (tt^A^)) C T z C°JT z X. 

Then J\(t) D J\(t) and J\(t) is a Lagrangian subspace of the symplectic 
space T Z C\/T Z \. In other words, J\(t) z = J\(t), where 

S z = {C G T Z C° A : a((,S) = 0}, S C T 2 . 

Given s G M \ {0}, sA is the singular extremal through sz G C\. Hence 
T z {Rz) C J\(t), Wt and Ja(^) C T z {Rz) z . This inclusion allows to make one 
more useful reduction. Namely, we set S z = T z (Wz) z /T z 'Rz, a symplectic 
space, dimE^ = 2(n — 3) and J\(t) is a Lagrangian subspace of S z . 
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Let L(E Z ) be the Lagrange Grassmannian: the manifold formed by all 
Lagrangian subspaces of S 2 . The curve t h- > J\(t) considered as a curve in 
L(E 2 ) is called the Jacobi curve associated to the extremal A. 

The dimension of a Lagrangian subspace is one half of the dimension of the 
ambient symplectic space. In particular, generic pair of Lagrangian subspaces 
has zero intersection. Jacobi curves are not at all generic, nevertheless, under 
very mild regularity assumption on the distribution (see [H]) they satisfy 
the following important property: J\{t) D J\(t) = for sufficiently small 
\t-r\^ 0. 

Let n tT be the linear projector of £ 2 on J\(t) along J\{t). In other words, 
T^tT '■ S 2 — > E z , 

7l tr\j x{t) =0, VTi T | jA(r) = 1. 

Lemma 1 (see|3j) We have: 

•where g\(t, r) is a symmetric function of (t, r) that is smooth in a neighbor- 
hood of (t, t), for all t out of a discrete subset of the domain of J\(-). 

In what follows, we tacitly assume that the value of t is taken out of the 
discrete subset provided by the lemma. A basic invariant of the parameter- 
ized singular extremal A is the generalized Ricci curvature 

t A (A(t)) d = f £ A (M). 

Generalized Ricci curvature depends on the parameterization of the ex- 
tremal; this dependence is controlled by the following chain rule. Let 
(p : K. — >• R be a change of the parameter, then: 

x Xov {\{<p{t))) = r A (A(^(t))) f{t) + (n - 3) 2 §(^), 

where §(</?) = — | ^f[^y^ > the Schwartzian derivative. The chain rule 
implies that the generalized Ricci curvature r A can be always made zero by 
a local reparameterization of the extremal A. We say that a local parameter 
t is a projective parameter if t\(t) = 0; such a parameter is defined up to a 
Mobius transformation. 



12 



Let f be a projective parameter, then the quantity: 

{dtf 

T=t 

is a well defined degree 4 differential on A; we called it the fundamental form 
on A. 

In arbitrary, not necessary projective parameter, the fundamental form 
has the following expression: 

Assume that A(X(t)) ^ 0, then the identity \A(X(s)) (^) | = 1 defines a 
unique (up to a translation) normal parameter s. 

Let z £ C\ and X s be the normally parameterized singular extremal 
through z. We set 

t(z) = t Xs (z), 

the projective generalized Ricci curvature. Then z h- > x(z) is a function on 
which depends only on A. 

Now come back to the case k = 2, n — 5. In this case, the fundamental 
form A is reduced to the famous Cartan's degree four form on the distribution 
constructed in [B] by the method of equivalence. The distribution is flat if 
and only if A = 0. 

Zelenko [14J performed detailed calculations for the rolling balls model. 
As before, let p be the ratio of the radii of the balls. We assume that 
1 < p < +oo. It appears that 

sign(v4) = sign(p — 3). 

Singular curves are just rolling motions along geodesies (i.e. along big cir- 
cles). The symmetry group acts transitively on the space of geodesies, hence 
function r must be constant in this case. We have: 

4y/35(p 2 + 1) 
r "3vV-9)(V-l)' 

In particular, the distributions corresponding to different p are mutually non 
equivalent and only the distribution corresponding to p = 3 is flat. 
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